WASP-37b: a 1.8 Mj exoplanet transiting a metal-poor star 
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ABSTRACT 

We report on the discovery of WASP-37b, a transiting hot Jupiter orbiting a m v = 12.7 
G2-type dwarf, with a period of 3.577469 ± 0.000011 d, transit epoch T = 2455338.6188 ± 
0.0006 (HJD^], and a transit duration 0.13041qqq} 8 d. The planetary companion has a mass 
M p = 1.80 ± 0.17 Mj and radius R p = L16+q% Rj, yielding a mean density of 1.15±^| pj. 
From a spectral analysis, we find the host star has M* = 0.925 ± 0.120 M Q , R* = 1.003 ± 
0.053 R , T cS = 5800 ± 150 K and [Fe/H] = -0.40 ± 0.12. WASP-37 is therefore one of the 
lowest metallicity stars to host a transiting planet. 

Subject headings: planetary systems — stars: individual: (WASP-37, GSC 00326-00658) - 
techniques: spectroscopic, photometric 



'Dates throughout the paper are given in Coordinated Universal Time (UTC) 
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1. Introduction 



Extrasolar planets show a huge diversity in their properties and this has important implications for 
theories of planet formation, structure and evolution. Systems with high orbital inclinations, in which the 
planet transits across the face of the host star as seen from Earth, are extremely valuable as they allow us to 
precisely measure many fu ndamental pla netary properties, including radius, mass and density, which can be 
used to test these theories (|Haswelll l2010). 



The parameter space which we are able to explore with transiting planets is biased by instrumental and 
observational limitations. However, many challenges faced by the current surveys are being overcome by the 
ability to decrease systematic noise and optimise follow-up strategies. Although the majority of the ~100 
transiting planets thus far discovered are short period, Jupiter-sized objects, they show a remarka ble variety 



in their physical and dynamical characteristics, such as the extreme eccentricity of HP 80606b (|N aef et al. 



20011 : iLaughlin et al.ll2009l : iMoutou et al 



2009 



Fossev et al.ll2009l : iGarcia-Melendo & McCullou gh 2009), 



the ultra-short period of WASP-19b (Heb b et all I201QT) and the puzzlingly low densities of WASP-17b 
jAnderson et allboid) and Kepler-7b (Latham et allboiol) . 



Here we describe the properties of a new transiting planet discovered by the SuperWASP survey, 
WASP-37b. SuperWASP has been a maj or contributor to the discovery of bright (9 < m v < 13) transiting 
planets since it began operation in 2004 (IPollacco et al.ll2006h . The project runs two stations, SuperWASP 
on La Palma, Canary Islands and WASP-S at SAAO in South Africa, each with a field of view of almost 
500 square degrees. A number of recent upgrades have been implemented to reduce systematic noise and 
improve photometric precision. These include reducing temperature fluctuations, which cause changes in 
the camera focus, by installing heating tubes, air conditioning and improving dome insulation. As a conse- 
quen ce, the varia tion in the stellar FWHM during the course of a night has been halved. For more details, 



see 



Barrosl feOIOh 



The planet host star WASP-37 resides in an equatorial region of the sky which is monitored by both 
WASP instruments, significantly increasing the amount of data collected on the target. It is accessible to 
observatories in both hemispheres and we present follow-up photometric and spectroscopic observations 
taken to establish the planetary nature of the transiting object and characterize it using the RISE (Liver- 
pool Telescope), Spectral (Faulkes Telescope South), SOPHIE (1.93-m OHP) and CORALIE (Swiss 1.2-m) 
instruments. 

This paper is structured as follows: Section [2] describes the observations, including the discovery data 
and photometric and spectroscopic follow-up. The results of the derived system parameters are presented in 
Section [3l including the stellar and planetary properties. Finally we discuss our findings in Section |4] 
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Fig. 1. — The combined, unbinned SuperWASP light curve for WASP-37, folded on the orbital period of P 
= 3.577 d. Superimposed is the model transit light curve, based on the system parameters determined from 
a global fit, see Section [3^21 

2. Observations 

2.1. SuperWASP photometry 

WASP-37 is a V = 12.7 mag star in the constellation Virgo, located at aj 2 ooo = 14 h 47 m 46 s .57, 
feooo = +01°03'53".9(GSC 00326-00658; 2MASS 14474655+0103538). It was observed by SuperWASP 
(La Palma) between March and June in 2008 and 2009, and by WASP-S (South Africa) during June to 
July 2008 and March to July 2009. A total of 22,593 photometric data points were obtained during these 
in tervals. The pipeline-processe d data were de-trended and searched for transits using the methods described 



in 



Collier Cameron et al.l (|2006T) . A periodic, transit-like signature was detected independently at the same 
period (3.58 d) in multiple cameras and in successive seasons, strongly suggesting that the transits were real, 
despite the faint magnitude. The folded light curve is shown in Figured] with the best-fit model. 

The star underwent several consistency tests aimed at eliminating false positives. It is isolated within 
the WASP photometric aperture in the Digitized Sky Survey (DSS) image , ruling out resolved blends. We 



also checked for signatures of unresolved blends (see lCameron et al.ll2007|) : No significant ellipsoidal vari 



ability was measurable in the WASP light curve folded on the transit period, nor any signs of a secondary 
eclipse. Moreover, the transit depth and duration yielded a planet-like radius for the companion, and a stel- 
lar density appropriate to a main-sequence host star of the effective temperature derived from the 2MASS 
colours. Having passed all these tests, the star was selected for follow-up observations. 



2.2. Photometric follow-up 



We obtained two further transit light curves of WASP-37 in order to refine the photometric parameters, 
and they are shown in Figure |2] All photometric data presented in this paper are available from the NStED 
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Fig. 2. — Photometry of transit events of WASP-37b on 2010 May 21 and 2010 June 30. The light curves 
have been offset from zero by arbitrary amounts for clarity. Superimposed are the model transit light curves 
based on the determined system parameters, see Section [3^21 Residuals from the fit are displayed underneath. 



database 

A partial transit was observed with LT/RISE on 2010 May 21. The full transit was not obtainable as 
the star set below the observing limits of the LT prior to egress. RISE is a frame transf er CCD located on 
the robotic 2 . 0-m L iverpool Telescope (LT) on La Palma with a broad band V + R filter ([Steele et al.l l2008: 
Gibson et al.ll2008r) . The telescope was defocused by —1.0 mm to give FWHM = 17 pixels = 9.2 arcsec, 
and the CCD was used in 2x2 binning mode with an exposure time of 65 s and effectively no dead time. 
This allowed 215 images to be taken over the 3.9 h p eriod, including 1 h o f observations before transit. The 
data were reduced using the ULTRACAM pipeline dDhillon et al. 2007 ) and differential photometry was 
performed relative to five nearby bright stars using an 18 pixel radius aperture. 



2 http://nsted.ipac. caltech.edu 
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A further partial transit was obtained on 2010 June 30 with the Spectrao camera (FS01) on the LCOGT 
2.0-m Faulkes Telescope South (FTS, Siding Spring, Australia). The Pan-STARRS-Z filter was used with 
the instrument in binning 2x2 mode, giving 0.303 arcsec/pixel, and no defocussing. The exposure time 
was 45 s, and 169 images were taken in the 3.5 h period. Data were reduced using a pipelineQ written 
at Liverpool John Moores University, then differential photometry performed using the IRAF/DAOPHOll 
package. Eight comparison stars were used with an 8 pixel radius aperture. Points with error bars larger than 
0.01 mag were removed from the second half of the dataset. These were likely caused by passing cloud, and 
the observation was stopped prematurely due to bad weather. 



2.3. Spectroscopic follow-up 



In order to establish the planetary nature and determine the stellar parameters, we obtained follow- 
up spectroscopic observations. Five spectra were taken between 2010 April 17 and 2010 May 28 with 
the stabilised echelle spectrograph SOPH IE at the 1 .93-m telescope of Observatoire de Haute-Provence 
( Perruchot et al. 2008 : Bouchy et al. 2009^ . The observations were all made with a signal-to-nois e ratio of 
S/N~20 in order to minimise the Charge Transfer Inefficiency (CTI) effect (|Bouchy et al.ll2009i) . Two 3 
arc-second diameter optical fibres were used, the first centred on the target and the second on the sky to 
simultaneously measure the background to remove contamination from scattered moonlight. A further 8 
spectra were obtained with the CORALIE Fibre-Fed Echelle Spectrograph on the Swiss 1.2-m telescope 
at ESO-La Silla, Chile between 2010 April 19 and 2010 June 29 with S/N~ 10-20 in dark/grey time to 
minimise contamina t ion from scatte red moonlight. The data were processed using the standard pipeline 
fearanne et al.lll996l : IPepe et alJEofCT) . 



The radial velocities (RVs) and line bisector spans (V^n, see lToner & Graylll988[) derived from cross 



correlation are shown in Table Q] and plotted with the best-fit model in Figure [3] No significant correlation 
is seen between the bisector span and radial velocity, with a supporting the sig nal's origin as a planetary 
companion rather than a blended eclipsing binary system (lOueloz et al] 12001). A Spearman rank-order 
correlation test also indicates that any correlation between the radial velocities and bisector is weak, with a 
probability of correlation of 0. 12. 



3. Results 



3 http://lcogt.net/en/network/2m 
4 http://telescope.livjm.ac.uk/Info/TelInst/Pipelines 

5 IRAF is distributed by National Optical Astronomy Observatories, operated by the Association of Universities for Research in 
Astronomy, Inc., under contract with the National Science Foundation, USA. 
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Table 1. Radial velocity (RV) and line bisector span (V^ pan ) measurements of WASP-37. 



BJD 


RV 


O-RV 


Vspan 


Instrument 


(2,400,000+) 


(kms- 1 ) 


(kms^ 1 ) 


(kms" 1 ) 




55304.4975 


7.892 


0.012 


0.024 


SOPHIE 


55305.4660 


8.223 


0.014 


-0.028 


SOPHIE 


55334.5449 


8.196 


0.019 


0.004 


SOPHIE 


55336.4764 


7.825 


0.021 


-0.088 


SOPHIE 


55345.4618 


8.115 


0.018 


0.013 


SOPHIE 


55305.7532 


8.152 


0.018 


-0.024 


CORALIE 


55310.8236 


7.678 


0.043 


-0.029 


CORALIE 


55311.8514 


7.919 


0.025 


-0.033 


CORALIE 


55312.8266 


8.064 


0.031 


-0.091 


CORALIE 


55321.7716 


7.668 


0.023 


0.058 


CORALIE 


55324.7114 


7.769 


0.022 


-0.011 


CORALIE 


55325.6868 


7.776 


0.024 


-0.025 


CORALIE 


55376.6684 


8.103 


0.031 


-0.004 


CORALIE 
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Fig. 3. — Upper panel: Phase folded radial velocity measurements of WASP-37, combining data from the 
Swiss 1.2-m/CORALIE (filled circles) and OHP/SOPHIE (open squares), and superimposed with the best-fit 
model RV curve based on the determined system parameters, see Section [3721 The centre-of-mass velocities, 
7sophie = 7.9714 kins" 1 and 7coralie = 7.9116 kms" 1 were subtracted from the RVs. The predicted 
Rossiter-McLaughlin effect is shown for a spin orbit aligned system with v sin i = 2.4 ±1.6 km s" 1 . Middle 
panel: Residuals from the orbital fit plotted against time. No long-term trend is visible. Lower panel: The 
bisector span measurements as a function of radial velocity, showing no correlation. The uncertainties in the 
bisectors were taken as twice the RV uncertainties. 
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Table 2. Stellar parameters of WASP-37 



Parameter (Unit) 



Value 



Photometric & spatial properties 
RA (J2000) 
DEC (J2000) 
V (mag) 
B - V (mag) 
J (mag) 
H (mag) 
K B (mag) 
/iRA (mas year -1 ) 
MDEC (mas year -1 ) 
[/(kms -1 ) 
^(kms- 1 ) 
^(kms" 1 ) 

Galactic longitude (/) (deg) 
Galactic latitude (b) (deg) 

Spectroscopic properties: 



Derived properties: 
A/* (M Q ) 
R* (R©) 
p* (pq) 
L* (L ) 
Age (Gyr) 
Distance (pc) 
Spectral type 



14:47:46.57 
+01:03:53.9 
12.704 ±0.149 
0.600 ± 0.050 
11.499 ±0.022 
11.181 ±0.026 
11.093 ±0.023 
-23.2 ±5.4 
22.8 ± 5.4 
-27.2 ±3.4 
7.8 ± 14.0 
45.5 ±3.4 
355 
52 



Teff (K) 


5800 ± 150 


logg (cgs) 


4.25 ± 0.15 


Vmic (kms -1 ) 


1.0 ±0.2 


"mac (kms -1 ) 


3.6 ± 0.3 


vsini (kms -1 ) 


2.4 ± 1.6 


[Fe/H] 


-0.40 ±0.12 


[Na/H] 


-0.38 ±0.16 


[Mg/H] 


-0.12 ±0.10 


[Si/H] 


-0.27 ±0.08 


[Ca/H] 


-0.25 ±0.13 


[Sc/H] 


-0.19 ±0.16 


[Ti/H] 


-0.28 ±0.14 


[V/H] 


-0.30 ±0.15 


[Cr/H] 


-0.37 ±0.16 


[Mn/H] 


-0.56 ±0.10 


[Co/H] 


-0.41 ±0.17 


[Ni/H] 


-0.37 ±0.08 


log A(Li) 


2.23 ± 0.13 


E(B - V) 


0.05 



0.925 ±0.120 
1.003 ± 0.053 

0.931^ 

0.953^ 
IV 



,+0.064 
-0.099 
3+0.184 
-0.114 
1+3 
-4 

343 ± 36 



G2V 



Note. — The photometric and spatial prop- 
erties are taken from or de rived using the fol - 
lowing sources: V (TASS, broege et all |2006|) . 
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B—V ICasagrande et alfcolch . J, H, K s (2MASS, 
Skrutskie et al. 2006), proper motions (NOMAD, 
Zacharias et alj 120050 . space velocities (see Sec- 
tion 13.1.11 and galactic co-ordinates (NED). The 
spectroscopic and derived properties are determined 
from the spectroscopic analysis described in Sec- 
tion [3~T] M* is the mean o f the masses f ound from 
the empirical relationship of lTorres et al] J2010h and 
stellar models; p* is found from the light curve ge- 
ometry; and the stellar radius is derived from p* and 
Af*. L* and age are determined from stellar mod- 
els. 
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3.1. Stellar parameters 



We used several techniques to infer the fundamental stellar properties of WAS P-37 and the resul ts are 
shown in Tabled We first performed a spectral analysis using the methods given in lGillon et all (120091) . The 
individual CORALIE and SOPHIE spectra from the standard reduction pipeline were co-added to produce 
a single spectrum with an average S/N of around 60:1. The H a line was used to determine the effective 
temperature (T c g), while the Na I D and Mg I b lines were used as surface gravity (log g) diagnostics. This 
yielded the values T e g = 580 ± 150 K a nd log g = 4.25 ± 0.15. The effective temperature indicates that 
the star is of spectral type G2 (lGrayll2008l) . 

The value for microturbulence (v m - lc ) was determined from Fe I using the iMagainl Jl984h m ethod. We 
assumed a value for macroturbulence (v mac ) of 3.6 ± 0.3 km s _1 based on the tabulation by Gray (2008), and 
an instrumental FWHM of 0.11 ± 0.01 A, determined from the telluric lines around 6300A. The projected 
stellar rotation velocity (v sin i) was determined by fitting the profiles of several unblended Fe I lines, giving 
a best-fitting value of v sin i = 2.4 ±1.6 km s -1 . 

Elemental abundances were determined from equivalent-width measurements of several clean and un- 
blended lines. The quoted errors include the uncertainties in T c g, log g and v m i c , as well as the scatter due 
to measurement and atomic data uncertainties. We obtain [Fe/H] = —0.4 ± 0.12 and log A(Li) = 2.23 ± 
0.13. Interstellar Na D lines are present in the spectra with equivalent widths of ^10 0-150 mA, indicating 
an extinction of E(B — V) = 0.05 using the calibrations of iMunari & Zwitten (119971) . 



Transit light curves provide a measure o f the mean stellar d ensity (p*) which can be used as a luminosity 
indicator for stellar evolutionary models (jsozzetti et al.l I2007I ). We used a Markov Chain Monte Carlo 
(MCMC) approach to globally model the photometric and radial velocity data (see Section I3T21 and obtained 
= 0.931^q'qqq p p). Va lues of T e g, [Fe/H] and p+ were compared with the theoretical stellar evolutionary 



models of lYi et all (120011) to obtain the following stellar properties: M* = 0.849^% M , R± = 0.977^°;°^ 



R©, L* = 0.9531q ^4 L , loggi so = 4.39^q'q3 and an age of 11^ Gyr. The model isochrones are shown in 
Figure |H 

We find that the measurement of logg from evolutionary models, loggi so = 4.39^Q'gg, is consistent 
with that found from the spectroscopic analysis, log g = 4.25 ±0.15. We re-fitted the isochrones using the 
lcr upper and lower limits of [Fe/H] and found that this affected the age and mass determinations somewhat, 
increasing the uncertainties in the age (8-13 Gyr) and stellar mass (0.806-0.914 M ). 

We computed the stellar distance using the following two methods. Firstly, the absolute bolometric 
magnitude, (calculated from the stellar luminosity) was compared to the extinction corrected V-band mag- 
nitude from TASS (m v = 12.7, Droege et all 20061). The bolo metric correction was taken as the range of 
values for a G2-type star given in ICoxl (bood) and bravl (I2008I) . This yielded a distance D = 338 ± 41 pc. 
Secondly, we compared the angular diameter O = 2 i?* / D = 0.025 ± 0.005 mas derived from the flux 
fitting method (see Section [3.1.1l > to the stellar radius, yielding a distance D = 360 ± 75 pc. It is reassuring 
that these values are similar, given that both methods rely on parameters from the stellar models and are 
therefore not independent. We take the distance to be the weighted mean of the two values, D = 343 ± 36 



j+0.02 
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Fig. 4. — Theoretical stellar isochrones from 



Yi et al 



1/3 

(2001) plotted as the density proxy R^/MJ against 



effective temperature for [Fe/H] = —0.4 and ages 1,9, 11, 14 and 16 Gyr (from left to right). The position 
of WASP-37 is indicated within the error ranges. 



pc. 



3.1.1. Independent checks 



The stellar mass can also be determined through an empirical calibration between T c g, p+ and [Fe/H] 
found from eclipsing binaries ( Torres et alJboiC : Enoch et aPboiO ). This yields a larger mass, M* = 1.002 



± 0.034 M than that from the evolutionary model fit, M* = 0.849+[J ^ M . A possible cause is that 
WASP-37 is on the edge of the calibration wher e there are no stars of similar T e g and [Fe/H] and log g. A 
similar situation was found for WASP-21, see lBouchy et all (120090 . In order to obtain a robust result, we 
have taken the stellar mass to be the average of the two values M*= 0.925 ± 0.120 M . 

We used archival magnitudes from 2MASS, DENIS, NOMAD, TASS and CMC 14 to derive an inde- 
pendent value of the stellar temperature using the flux-fitting method ([Morossi & Malagninilll985h . We find 
leffflux = 5840 ± 140 K, in excellent agreement with the temperature derived from spectral analysis. 

Independent measures of v sin i and [Fe/H] can be derived from the sh ape of the spectral cross correla- 
tion function (CCF). Using the empirical calibrations of iBoisse et all (120101) for the SOPHIE HE mode, and 
inputing Contrast = 4.4 ± 1. 3% and FWHM = 9.48 ± 0.06 kms" 1 from the SOPHIE CCF and B — V = 
0.60 ± 0.05 (calculated using ICasagrande et al hoick we find [Fe/H]ccF = —0.32 ± 0.17 and v siniccF = 
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3.2 ± 1.0 km s . These are also in excellent agreement with the values found from spectral analysis. 



The stellar evolution modelling provides colour indices which can be compared to observational values. 
We corrected the 2MASS magnitudes for extinction using A(V ) = 3.1 xE(B — V), appropriately scaled 
for each passband following ICardelli. Clayton. & Mathisl (119891) . The model colours were transformed from 
the ESO to the 2MASS system following Carpenter "(feoOll ). The 2MASS and model colours agree well, as 
shown in Table [3] 

The large relative error on v sin % provides little constraint on t he maximum stellar rotation period 
P rot < 2irR*/v s'mi < 21^g 4 d, or gyrochronological age; 1-28 Gyr JBarnesll2007h . Equally, the lithium 



abundance only implies a lower limit on the age of 1-2 Gyr (|Sestito & Randichl l2005). 



Given the low metallicity of WASP-37, we inve stigated whether it could be a member of the thick 
disk using the method outlin ed in iBouchy et all (I2010h . We calculated the space velocities according to 
Johnson & Soderbloml (119871) . using the proper motions, distance and 7 given i n Tables [2] and [H T he ve - 
locities were updated to J2000 and corrected for solar motion using data from iDehnen & Binneyl (11998b 
giving U = -27.2 ± 3.4 km s -1 , V = 7.8 ± 14.0 km s" 1 and W = 45.5 ± 4.4 km s _1 (right-hand coordinate 
system). 

In order to quantify the probability that WASP-37 is a member of the thick disk, we c ompared its 



metal licity and space velocities to the properties of similar stars in the Besancon Galactic model (IRobin et al. 



2003). We generated 78,000 model stars in 20 realisations and returned all stars with apparent magnitudes 
between m v = 8-16, spectral type F5-K9 and luminosity class IV and V for 30 square degrees in the 
direction of WASP-37 (galactic coordinates I = 356°, b = 53°, NED). For the stars with [Fe/H] and space 
velocities within 3a of WASP-37, half were members of the old thin disk (>6 Gyr), one third were part 
of the thick disk, and the remainder were 1.5-6 Gyr thin disk members. Abundance analysis shows that 
WASP-37 has enhanced alpha elements and although this is not accounted for in the simulation, we estimate 
that it increases the probability of thick disk membership but that it is also equally likely to be part of the old 
thin disk. It is reassuring that this independent analysis also finds that WASP-37 is part of an old population. 



Table 3. Derived colors of WASP-37 





J-K 


J - H 


H-K 


Model 


0.377 ±0.026 


0.293 ± 0.026 


0.019 ±0.010 


2MASS 


0.380 ±0.032 


0.304 ± 0.034 


0.076 ±0.035 



Note. — The model colors are derived from stellar evolutionary 
models and transformed to the 2MASS system. The 2MASS colours 
are corrected for extinction. See Section l3.1.1l for details. 
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Table 4. System parameters of WASP-37b 



Parameter (Unit) 



Value 



Photometric parameters: 

P(d) 
T (HJD) 

Tdur (d) 

AF = (_Rp/R») 2 

Rp/R* 

a/R* 

b 

Spectroscopic parameters: 
if (kms^ 1 ) 
7SOPHIE (kms^ 1 ) 
7CORALIB (kms -1 ) 
e 

Derived parameters: 
M p (Mj) 

Hp (Rj) 
pp(pj) 

log g P (cgs) 
a(AU) 
2eq,A=0 (K) 



3.577469 ±0.000011 

2455338.6188 ± 0.0006 

1304+ - 0018 
u.iout_ 0017 

01427+ ' 00030 

11946+ 00126 
U.115"K)_ 00096 

9.567 ± 0.648 



0.198i°;ig 

,,+0.77 
-0.86 



88.82] 



0.2507 ± 0.0084 
7.9714 ± 0.0024 
7.9116 ±0.0024 
(fixed) 



1.80 ±0.17 
1.161 
1.15j 



+0.07 
0.06 
+0.12 
0.15 
+0.03 



3.48_ 04 
0.0446 ± 0.0019 

,+25 
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Table 5. Limb darkening coefficients 



Light curve 


Band 


ai 


CL2 


a 3 


0,4 


SuperWASP 


R 


0.643 


-0.278 


0.800 


-0.413 


RISE 


V 


0.570 


-0.131 


0.791 


-0.420 


FTN 


Z 


0.623 


-0.331 


0.672 


-0.346 
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3.2. Planet parameters 



To determine the properties of the planet, we simultaneously modelled the SuperWASP, LT and FTN 
light curves and the C ORALIE and SOPHIE rad ial velocities with a global MCMC fit. Details of this 
process are described in lCameron et all (12007b and lPollacco et all (120081) . The free parameters in the fit are: 



orbital period P; transit epoch To; transit duration 7d ur ; ratio of planet to star radius (R p /R*) ; impact 
parameter b; RV semi-amplitude K ; Lagrangian elements e cos ui and e sin u where e is the eccentricity and 
uj is the longitude of periastron; and the systematic offset velocity 7. In this particular case, two systematic 
velocities were fit to allow for instrumental offsets between the SOPHIE and CORALIE datasets. 

We propagated the uncertainty in the stellar mass through to the deriv ed parameters (M* = 0.925 ± 
0.120 M , see Section [3.1.1l >. We used the 4 coefficient Iciaret ( 200C , 20041) non-linear limb darkening co- 
efficients for the appropriate stellar temperature and photometric passband appropriate for each light curve, 
see Table [5] The limb darkening coefficients are re-calculated at each step in the MCMC chain to take into 
account the uncertainties in the stellar temperature and radius. 

The photometric error bars were scaled to account for any underestimations so that the best fitting 
model results in x 2 eA = x 2 /dof = 1 (dof = number of points — number of fitted parameters). The radial 
velocity error bars did not require rescaling with a jitter term in order to obtain x 2 c d = 1- The RMS of the 
residuals from the best-fit model are: RMS wasp = 0.0292 mag, RMS L t = 0.0013 mag, RMSfts = 0.0025 
mag, RMSsophie = 0.0194 kms _1 and RMScoralie = 0.0380 kms -1 . By removing the outlying point at 
8.064 km s _1 , we obtain RMScoralie = 0.0164 km s -1 . 

An initial fit was performe d in which the eccentri city was allowed to float, yielding e = 0.036 ± 0.021. 
From the Lucy-Sweeney test (ILucy & Sweeneyill97ll) . we find it is only significant at the la level, hence 



a circular orbit was adopted (x 2 i rc = 18.3, x 2 cc = 14.0 with 13 RVs and 5 free parameters: K, 7coralie, 
7SOPHIE, e cos u and e sin u). We determined an upper limit for the eccentricity as e upP er = e+3a = 0.078. 
We also fitted a linear trend in the RVs to search for a third body in the system, but found that the decrease in 
X 2 did not warrant the extra free parameter. From the F-test, we found the significance of the trend was less 
than la (x 2 = 18.3, x 2 ren d = 15.4 with 13 RVs and 4 free parameters: K, 7coralie, 7sophie and linear 
trend 7). 

We performed the final MCMC fit with a chain length of 20,000 points and the resulting best-fit pa- 
rameters and uncertainties are shown in Tableg] WASP-37b has P = 3.577469 ± 0.000011 d, M p = 1.80 ± 



0.17 Mj, R p = 1.16l°;°6 R J and Pp = l - 15 -o.lt PJ- 



+0.12 



4. Discussion 



We report the discovery of a new transiting hot Jupiter, WASP-37b. The planet has a radius comparable 
to Jupiter (.Bp = 1.16 + ££j Rj) but is more massive (M p = 1.80 ± 0.17 Mj) and hence is very similar to 



WASP-5 (And erson et al. 



2008: ISouthworth et al.ll2009h . The surface gravity of WASP-37b is high com- 



pared to other planets with similar orbital periods (log g p = 3.48 



+0.03 



^2.0 



-0.04' yp - -' u - z -2.8 



ms 



K P P = 1.151 



0.12 

0.15 
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Fig. 5. — Evolutionary mass tracks from lYi et al.l (|2001h for [Fe/H] = —0.4 and stellar masses of 0.95, 0.9, 
0.85 and 0.8 M (from left to right). The current position of WASP-37 is shown as a square within the la 
uncertainty ellipse. The 12, 13, 14, 15 and 16 Gyr points are shown on the best-fit, 0.85 M track. WASP-37 
is predicted to evolve onto the subgiant branch at age ~15 Gyr and red giant branch ~1 Gyr later. 



pj), which lies significantly above the period-gravity correlation proposed by lSouthworthl (I2010f) . 



By comparison with the theoretical models of Fortney. Marley. & Barnes! (12007b and lBaraffe. Chabrier. & Barman 



(|2008h . we find that WASP-37b has an inflated radius. The models do not cover the age range of WASP-37, 
however given that planetary radii are thought to decrease with age, the expected radius of WASP-37b will 
be smaller than those given for the oldest model. For an equivalent semi-major axis with solar insolation of 
0.445 AU, the measured mass and radius do not suggest the presence of a core of heavy elemen ts for models 



older than 1 Gyr, and this i s consistent with the correlation between core mass and metallicity (iGuillot et al. 
2006: [Burrows et al.ll2007r) given the low stellar metallicity. 



We detect no significant orbital eccentricity and find an upper limit of 0.098. Nor do we find evidence of 
any long-term trend caused by a third body. Our radial velocity observations span ~ 70 days and continued 
monitoring would reveal the presence of any longer-term trends and further constrain the eccentricity. 

With a stellar effective temperature of 5800 ± 150 K and orbital period of 3.577469 ± 0.000011 d, 
the planet is predicted to have a equilibrium temperature of 1323l^| K (assuming zero albedo and full 
redistribution of the incident radiation). Secondary eclipse measurements of WASP-37 may be a way of 
constraining the influence of metallicity on the presence or not of atmospheric temperature inversions. 



The relatively large uncertainties in the stellar parameters are caused by the low S/N of the combined 
spectrum used for spectral analysis, due to the faintness of the target (m v = 12.7). Further spectra would 
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Fig. 6. — Histogram of the metallicity of transiting planet host stars and stars w ith and without detections of 
RV planets from the HARPS GTO "high-precision" sample dsousa et allboOof ). WASP-37 is clearly in the 
tail of the metallicity distribution for stars hosting planets. 



help to constrain the spectroscopic and derived parameters. A more precise value of v sin i could also 
be found from an observation of the Rossiter-McLaughlin (RM) Effect (IRossiten 1 1 924l : iMcLaughlinl 1 1 924l : 



Gaudi & Winn 



2007), which could provide an independent constraint of the age through gyrochronological 



relations ( Barnesll2007 ). The slow rotation (v sini = 2.4 ±1.6 kms x ) and faintness of the host star would 
make observations of the RM effect challenging but possible given the predicted amplitude of 35 ms -1 . 

We find a small discrepancy between the stellar mass found from evolutionary models and an empirical 
calibrations based on eclipsing binaries, and take a mean of the two values in this analysis. WASP-37 is less 
massive than the Sun, M* = 0.925 ± 0. 120 M , and is one of the older stars hosting a transiting planet at ~1 1 
Gyr. Although the age is relatively uncertain, it suggests that giant planet formation was taking place when 
the Milky Way was still relatively young. The model isochrone, shown in Figure [51 suggests that WASP-37 
is due to evolve from a main-sequence star into a subgiant at an age of ~15 Gyr and subsequently move on to 
the giant branch after a further ~1 Gyr. Although no meaningful constraint can be placed on t he stellar tidal 
qualit y factor, Q s , we note that if is greater than 10' then the planetary spiral-in timescale qLevrard et al. 
2009) is longer than the main sequence lifetime, and the planet will be engulfed by the expansion of the 
stellar radius during the giant phase less than 0.5 Gyr after the star has evolved onto the giant branch. 

WASP-37 ha s the second lowes t metallicity ([Fe/H] = —0.4 ± 0.12) of the transiting planet host stars 
after WASP-21 (|Bouchy et al.ll2010l) . and is firmly in the tail o f the metallicity distribution fo r exoplanet 
host stars as shown in Figure [6] The figure, similar to that in lAmmler-von Eiff et al.1 (|2009l ). shows the 
metallicity distributions for stars with transiting planets and stars with and without detections of RV planets 



-17- 



from the HARPS GTO "high-precision" sample (iSousa et al.ll2008f) . The stars hosting transiting and RV 
planets show a very similar distribution, whereas the non-planet hosts appear to extend to m uch lower 
metallicities, suggesting that the heavy metal con tent plays a significant role in planet formation (Gonzalez 



1998; Santos et al. 2004; Fischer & Valenti 2005). 



The SuperWASP Consortium consists of astronomers primarily from Queens University Belfast, St 
Andrews, Keele, Leicester, The Open University, Isaac Newton Group La Palma and Instituto de Astrofsica 
de Canarias. The SuperWASP-N camera is hosted by the Issac Newton Group on La Palma and WASP- 
South is hosted by SAAO. We are grateful for their support and assistance. Funding for WASP comes from 
consortium universities and from the UK's Science and Technology Facilities Council. FPK is grateful to 
AWE Aldermaston for the award of a William Penny Fellowship. Based on observations made at Observa- 
toire de Haute Provence (CNRS), France and at the ESO La Silla Observatory (Chile) with the CORALIE 
Echelle spectrograph mounted on the Swiss telescope. The RISE instrument mounted in the Liverpool Tele- 
scope was designed and built with resources made available from Queens University Belfast, Liverpool John 
Moores University and the University of Manchester. The Liverpool Telescope is operated on the island of 
La Palma by Liverpool John Moores University in the Spanish Observatorio del Roque de los Muchachos 
of the Instituto de Astrofisica de Canarias with financial support from the UK Science and Technology Fa- 
cilities Council. The research leading to these results has received funding from the European Community's 
Seventh Framework Programme (FP7/2007-2013) under grant agreement number RG226604 (OPTICON). 
We thank Tom Marsh for the use of the ULTRACAM pipeline. 
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